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I.  INTRODUCTION 


The  demand  for  optimum  maneuverability  in  the  design  of 
tactical  missiles,  as  well  as  spacecraft  launch  vehicles, 
necessitates  the  application  of  active  control  systems.  The 
control  forces  and  momentum  generated  by  such  systems  are 
essential  to  accomplish  specified  trajectory  changes. 

Of  the  various  methods  of  active  trajectory 
correction,  thrust  vector  control  (TVC)  systems  offer  the 
only  means  of  trajectory  control  that  are  independent  of  the 
external  forces  on  the  vehicle.  Such  capability  is  required 
when  the  flow  past  the  vehicle's  external  aerodynamic 
surfaces  (canards,  wings,  tails,  etc.)  is  insufficient  to 
generate  the  necessary  control  forces.  This  commonly  occurs 
during  low-speed  flight,  such  as  at  launch  or  during 
hovering.  High  angle  of  attack  flight  may  also  lead  to 
regimes  where  conventional  aerodynamic  lifting  surfaces  are 
inadequate  [Ref.  1].  Furthermore,  some  design 
configurations  render  external  control  surfaces  impractical, 
such  as  tube-launched  missiles. 

Three  different  methods  of  TVC  have  been  developed  and 
operationally  tested;  movable  no?.2le,  secondary  injection 
(internal  fluid  injection),  and  mechanical  jet 
deflection  systems  [Refs.  2,3).  Comparison  studies  favor 


the  latter  category  which  are  characterized  by  relatively 
low  actuation  torques.,  rapid  response  and,  in  the  case  of  a 
jet  vane  system,  a  small  installation  envelope.  Of  the 
various  TVC  systems,  jet  vanes  have  proven  to  be  the  most 
effective  in  delivering  large  thrust  vector  deflection 
angles  (up  to  30  degrees.)  [Ref.  3].  They  also  offer  the 
possibility  of  roll  control  and,  for  the  case  of  a  single 
exhaust  nozzle  arrangement,  a  jet  vane  device  is  the  only 
practical  system  which  could  perform  the  task. 

There  are,  however,  disadvantages  in  the  use  of  jet 
vane  TVC.  First,  there  are  thrust  losses  on  the  order  of 
3-5%  with  the  use  of  non-retractable  vanes.  In  addition, 
the  ability  to  achieve  such  high  thrust  deflection  angles 
can  result  in  axial  trust  losses  proportional  to  50  -  100% 
of  the  resulting  side-force  (Ref.  3].  Finally,  the  most 
significant  problem  with  the  use  of  jet  vanes  is  the  large 
thermal  loading  they  experience  during  exposure  to  high 
temperature,  supersonic,  particulate-laden  rocket  exhaust 
gases.  In  the  past,  this  problem  has  limited  jet  vanes  to 
short  duration  use  in  engines  with  low  temperature  non- 
metal  ized  propellants. 

Jet  vane  TVC  dates  back  more  than  50  years  to  the 
rockets  designed  by  R.H.  Goddard,  and  has  extended  to 
Redstone,  Sargent,  Tales,  Pershing,  and  Anglo  II  and  III 
engines,  as  well  as  several  installations  in  smaller 
tactical  rockets  (Refs.  1,2].  Due  to  the  material  problems 
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associated  with  the  use  of  high  energy  aluminized  propel¬ 
lents  and  the  advance  of  other  TVC  systems,  jet  vane  devices 
had  been  furthered  studied  only  on  a  low  level.  New  mission 
requirements,  however,  coupled  with  the  above  mentioned 
system  capabilities  have  resulted  in  a  renewed  level  of 
interest  in  the  application  of  jet  vane  TVC  especially  for 
launch  phase  trajectory  control. 

The  aerodynamic  design  of  the  jet  vane  can  be  performed 
with  reasonable  accuracy  on  the  basis  of  supersonic  flow 
theory  with  boundary  layer  corrections  [Pefs.  3,4].  It  is 
significantly  more  difficult  to  calculate,  with  the  same 
degree  of  certainty,  the  vane  heat  transfer  characteristics 
due  to  its  severe  thermal  environment.  This  inability  to 
accurately  predict  jet  vane  transient  thermal  response 
results  in  design  material  selection  and  configuration  based 
solely  upon  past  experience  and  costly  experimental 
findings. 

Hence  the  development  of  a  dynamic  computational  model 
capable  of  calculating,  with  the  requisite  design  certainty, 
the  transient  thermal  response  of  such  systems  would  truly 
be  beneficial.  Such  a  model  could  provide  a  capacity  for 
design  optimization  which  is  virtually  non-existent  in 
current  design  efforts. 


II.  BACKGROUND 


The  basic  goal  of  jet  vane  heat  transfer  studies  is  to 
develop  a  capacity  to  accurately  predict  the  transient 
thermal  behavior  of  the  vane.  The  severe  thermal 
environment  in  which  the  vane  operates  makes  analysis  and 
modeling  of  the  energy  transport  processes  extremely 
difficult.  A  complete  analysis  would  have  to  account  for 
multi-phase,  multi-component,  three-dimensional  and  time 
dependent  effects  in  the  presence  of  shocks,  boundary 
layer  transition,  and  turbulence,  separated  flows,  surface 
ablation,  chemical  reaction,  solid-body  and  gaseous 
radiation  [Ref.  5], 

It  can  be  assumed  that  all  the  details  of  the  energy 
transport  process  will  never  be  fully  defined.  The  net 
effect  of  all  the  various  complications  mentioned  above  is 
to  transfer  energy  to  and  from  the  vane  system.  This 
process  causes  energy  flow  within  the  vane  which  in  turn 
results  in  a  vane  temperature  distribution  that  varies 
with  time. 

This  transient  thermal  response  should  reflect  the 
nature  of  the  surface  heat  transfer  processes  that  drive 
it,  as  well  as  the  thermal  impact  of  various  vane  design 
parameters.  Hence  the  transient  response  is  the 


4 


"signature"  of  the  combined  effects  of  both  the  thermal 
environment  and  the  vane  construction  [Ref.  5].  This 
concept  is  the  basis  for  the  modeling  and  simulation  (M&S) 
study  described  herein. 

The  goal  of  the  M&S  study  has  been  to  develop  a 
dynamic  model  of  vane  thermal  behavior.  A  model  that 
could  predict,  with  sufficient  accuracy,  the  transient 
thermal  response  of  the  vane  would  be  an  invaluable  design 
tool.  For  instance  it  could  be  used  to  predict  the 
effects  of  design  changes  without  the  need  for  costly 
experimentation  and  testing.  Also,  sub-scale  test  data 
could  be  used  to  predict  the  full  scale  jet  vane  response. 
Most  significant  is  the  model's  inherent  ability  to  easily 
accommodate  design  optimization.  The  model  could  be  used 
for  deductive  purposes  also.  For  example,  measured  vane 
temperature  data  could  be  used  for  parameter  estimation. 
That  is,  use  of  the  model  to  deduce  what  would  have  had  to 
occur  in  order  to  obtain  the  measured  results. 

Furthermore,  estimation  of  local  temperatures  could  be 
made  at  points  in  the  vane  where  accessibility,  sensor 
survivability  and  other  conditions  render  measurement 
impractical  if  not  impossible. 

Research  has  shown  that  the  traditional  approach  to 
modeling  a  system  such  as  this  would  be  to  construct  a 
comprehensive  model  that  would  treat  the  flow  environment 

of  the  vane  and  the  vane  itself  in  fine  numerical  detail. 
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Numerical  techniques  such  as  finite-element  analysis  or 
the  finite-difference  method  would  then  be  employed  to 
determine  both  the  steady-state  heat  transfer  and  the 
temperature  distribution.  The  M&S  method  can  thus  be 
considered  a  vast  simplification  of  a  conventional 
numerical  model.  A  basic  assumption  of  this  method  is 
that  sufficient  accuracy  can  be  obtained  if  the  flow  and 
vane  are  considered  to  be  made  up  of  relatively  few 
thermal  components.  An  important  facet  of  the  M&S  model 
is  the  computation  of  temperature  distribution  as  a 
function  of  time.  This  transient  response  capability 
lends  itself  well  to  model  verification  based  on  actual 
test  data. 

The  modeling  procedures  used  in  this  study  are  based 
on  the  concept  that  the  thermal  processes  associated  with 
jet  vane  operation  can  in  fact  be  put  into  a  transfer- 
function  form  commonly  found  in  the  discipline  of 
automatic  controls.  Hence,  experimental  data  can  be 
related  to  the  various  gains  and  frequencies  of  the 
transfer  function  by  well  established  system  identifi¬ 
cation  methods.  The  establishment  and  verification  of 
these  relationships  has  been  the  objective  of  this  study. 
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III.  THERMAL  MODELING  APPROACH 


A.  CONCEPT 

The  modeling  procedures  used  in  this  study,  as  pre¬ 
viously  mentioned,  are  based  on  the  concept  that  the 
transient  thermal  response  of  a  jet  vane  can  be  formulated 
by  the  use  of  transfer  functions.  These  transfer  functions 
are  developed  from  equations  defining  the  various  heat 
transfer  processes  occurring  during  jet  vane  operation. 

The  "  lumped-heat-capacity '•  method  is  the  analytical 
approach  employed  to  develop  such  equations.  This  analysis 
yields  reasonably  good  results  when  the  resistance  to  heat 
transfer  by  conduction  is  small  compared  with  the 
convective  resistance  at  the  surface.  Hence,  the  non- 
dimensional  Biot  number  (Bi)  becomes  a  most  significant 
system  parameter.  The  Biot  number  denotes  the  resistance 
to  conduction  in  the  vane  relative  to  the  convective 
resistance  of  the  rocket  exhaust  flow.  Therefore,  to 
ensure  that  the  error  resulting  a  lumped  capacitance 
analysis  is  minimized,  the  following  condition  must  be 
satisfied, 

Bi  ®  hV/KA  <0.1  (3.1) 

where  h  is  the  heat  transfer  coefficient,  K  is  the  thermal 
conductivity  and  V  and  A  are  the  volume  and  surface  area 
respectively  [Ref.  6;pp.  134-135]. 
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B .  METHODOLOGY 


A  jet  vane  of  a  particular  design  was  hypothetically 
configured  to  consist  of  several  discrete  thermal  parts  or 
"lumps'*.  Then  a  thermal  circuit  was  constructed  with  nodes 
located  at  the  center  of  mass  of  each  lump.  Figure  3.1  is 
a  simplified,  two  lump,  example  of  such  a  circuit. 

As  evidenced  by  Figure  3.1,  thermal  networks  consist  of 
only  two  types  of  components;  thermal  resistances  (R)  and 
thermal  capacitances  (C) .  The  type  of  thermal  resistance 
modeled  depends  on  the  mode  of  heat  transfer  being 
considered.  For  conduction  between  nodes, 

R  =  L/KA  (3.2) 

whereas  for  convection  to  the  vane, 

R  **  1/hA  (3.3) 

Similarly,  there  are  resistances  associated  with  radiation 
which  is  a  subject  treated  in  subsequent  sections. 

Thermal  capacitance,  on  the  other  hand,  is  independent 
of  the  mode  of  heat  transfer,  and  is  defined  as, 

C  =»  pVca  (3.4) 

where  p  is  the  density  of  the  vane  material  and  cs  is  the 
associated  specific  heat. 
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With  the  thermal  network  thus  defined,  the  requisite 
governing  equations  are  developed  by  means  of  an  energy 
balance  at  each  node.  An  energy  balance  is  defined  as 
follows: 

Rate  of  heat  Rate  of  heat  Rate  of  energy 

flow  into  =  flow  out  of  +  storage  at  a 
node .  node .  node . 

For  instance,  a  balance  at  node  1,  of  Figure  3.1, 
would  yield: 


T 


r 


R1 


clsTl 


(3.5) 


where  the  letter  s  is  the  Laplace  variable. 

The  nodal  temperature  can  be  expressed  explicitly  as 
follows: 


T 


1 


'n 


1  +  rS 


(3.6) 


where, 


111 

R„  Rx  R2  (3-7) 

and  the  nodal  time  constant  is  given  by 

t  *  RnC1  (3.8) 

A  heat  balance  at  node  2  would  yield  similar  results, 
thereby  providing  sufficient  information  to  determine  both 
nodal  temperatures  as  a  function  of  time. 
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The  computational  method  and  the  results  of  the 
application  of  the  modeling  approach  outlined  above  are 
described  in  later  sections. 

C .  APPLICATION 

The  modeling  approach  described  above  can  be  applied 
to  various  TVC  jet  vane  configurations.  In  this  study, 
however,  a  single  design  was  used  as  the  basis  for  all 
models.  The  particular  design  chosen  was  the  retractable 
jet  vane  employed  in  the  Stowable  Three-Axis  Reaction 
Steering  (STARS)  System,  Figure  3.2,  currently  under 
development  at  the  Naval  Weapons  Center  (NWC) ,  China  Lake 
[Ref.  7] 

The  STARS  jet  vane  is  constructed  of  copper 
impregnated  tungsten  (10%Cu/W),  comprised  of  10%  copper  and 
90%  tungsten  (by  volume) ,  and  the  vane  support  hardware  is 
steel.  A  major  operational  feature  of  the  STARS  system  is 
the  variable  positioning  of  the  TVC  vanes  relative  to  the 
radial  plane  of  the  rocket  nozzle  exit.  The  TVC  vane 
control  system  currently  being  tested  is  capable  of  moving 
the  vanes  from  the  fully  inserted  "active"  position,  to  the 
fully  retracted  "stowed"  position  as  illustrated  in  Figure 
3.3.  It  should  be  noted,  however,  that  all  the  analyses 
performed  in  this  report  were  based  solely  on  the  TVC  vane 
remaining  in  the  "active"  position. 
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TVC  VANE 


WITH  DCOICATEO  ACTUATOR  ANO  PNEUMATIC 
RETRACTION  DEVICE 


Figure  3.4  shows  a  single,  full-scale  TVC  jet  vane 
mounted  on  a  mock-up  of  a  Vertical  Launch  Modular  Booster 
(VLMB)  rocket  motor  and  nozzle  assembly.  Figures  3.5  and 
3.6  are  used  to  demonstrate  the  TVC  vane  in  the  actual 
"stowed"  and  "active"  positions  respectively. 


IV.  MODEL  DEVELOPMENT 


A.  BASIC  DESCRIPTION 

Investigations  into  the  feasibility  of  applying  the  M&S 
approach  to  determine  the  dynamic  thermal  characteristics  of 
the  STARS  jet  vane  system  began  in  early  1986  [Ref.  5].  The 
results  of  these  preliminary  studies,  as  outlined  below, 
form  the  basis  of  the  work  reported  earlier.  The  thermal 
model  described  in  this  chapter  is  but  one  product  of  these 
early  investigations. 

The  current  model  considers  only  two  thermal  energy 
input  processes:  forced  convection  at  the  vane  surface  and 
stagnation-point  heat  transfer  at  the  leading  edge.  These 
processes  are  driven  by  the  flow  stagnation  temperature 
which  is  a  function  of  the  rocket  motor  thrust  levels.  The 
basic  assumption  is  that  the  stagnation  temperature  will 
follow  the  rise  and  fall  of  the  thrust  without  a  significant 
dynamic  offset.  Hence,  the  basic  output  of  the  model  is  a 
series  of  functions  that  give  local  vane  temperatures  as  a 
function  of  thrust. 

several  discretization  schemes  were  tested  and  it  was 
found  that  four  H lumps”  were  sufficient  to  indicate  the 
general  thermal  behavior  of  the  vain.  The  four  lumps  are: 
the  vai  a  tip,  vane  body,  shaft  and  mount. 
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B.  GEOMETRICAL  AND  THERMAL  ESTIMATIONS 

In  one  step  of  the  preliminary  study-  the  actual  vane 
design  was  configured  as  a  collection  of  lumps  with  regular 
geometries  more  suitable  for  calculating  the  thermal  con¬ 
duction  properties.  Figure  4.1  illustrates  the  actual  vane 
geometry  design  while  the  discretized  vane  is  shown  in 
Figure  4.2  [Ref.  5].  The  fictitious  thermal  vane  was 
hypothesized  as  consisting  of  three  rectangular  solids: 
tip,  fin  and  shaft.  The  tip  was  arbitrarily  sized  so  as  to 
have  a  chord  length  of  1.0%  of  the  total  vane  chord  of  3.75 
inches.  With  the  length  fixed,  the  tip  was  shaped  so  as  to 
have  the  same  lateral  area  (chord  x  thickness)  as  that  of 
the  actual  geometry.  The  discretized  tip  was  separately 
identified  to  account  for  the  stagnation  properties  of  the 
convection  heat  transfer  near  the  vane  leading  edge. 

The  remaining  portion  of  the  discretized  fin  was  assumed 
to  be  subjected  to  thermal  convection  of  a  turbulent 
boundary  layer  type.  This  section  was  lumped  into  a 
rectangular  solid  of  thickness  equal  to  the  average  value  of 
the  tapered  fin  minus  the  tip  and  a  span  equal  to  the  span 
of  the  actual  fin. 

The  remaining  length  dimension  was  established  by 
equating  the  volume  of  the  rectangular  solid  to  that  of  the 
actual  fin  (with  tip  removed) .  The  vane  shaft  was  similarly 
modeled  into  a  rectangular  solid,  and  assumed  to  be 
subjected  to  conduction  heat  transfer  only. 
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MOUNT 


Figurtt  4.2.  Diacratiaad  Version  of  TVC  Jot  Vana 
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Tne  vane  support  mount,  as  shown  in  figure  3.2,  is  a 
relatively  complex  structure.  In  order  to  pursue  the  basic 
feasibility  of  the  M&S  method  [Ref.  5] ,  no  attempt  was  made 
to  model  the  thermal  resistance  of  this  component  in  detail. 
Instead,  the  conduction  length  and  cross-section  for  the 
main  mount  components  were  estimated  using  a  rough  scalling 
procedure.  An  appropriate  thermal  resistance  for  the  mount 
was  then  computed  from  an  analysis  of  the  analogous 
electrical  circuit.  In  the  preliminary  studies  the  mount 
thermal  resistance  was  used  as  an  “adjustable  parameter"  in 
seeking  agreement  with  the  available  experimental  data. 
Hence,  the  determination  of  the  correct  value  for  the  mount 
thermal  resistance  is  one  of  the  goals  of  the  system 
identification  procedure  described  in  later  chapters. 

A  thermal  node  was  located  at  the  center  of  mass  of 
each  of  the  components  described  above.  These  nodes  are 
assumed  to  be  the  location  of  the  energy  storage  associated 
with  the  entire  mass.  The  thermal  capacity  of  the  mount 
was  considered  to  be  infinite,  thus  allowing  its  nodal 
temperature  to  remain  constant  at  ambient  temperature 
(ground)  during  the  simulation. 

It  should  be  noted  that  the  primary  goal  of  the 
discretization  process  was  to  obtain  adequate  agreement  with 
test  results  using  the  minimum  number  of  thermal  components. 
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This  logic  leads  to  the  necessary  simplicity  that  is  key  to 
the  M&S  approach  and  the  results  to  be  described  later 
justify  the  development  of  such  "simple"  models. 


C.  DEVELOPMENT  OF  THE  GOVERNING  EQUATIONS 

Applying  the  methodology  described  in  the  previous 
chapter  to  the  discretization  process  outlined  above  yields 
thermal  model  illustrated  in  Figure  4.3.  The  nominal  values 
of  the  various  thermal  components  listed  in  the  figure,  were 
developed  in  the  preliminary  studies  of  Reference  [8].  The 
governing  equations  can  thus  be  formulated  by  using  the 
energy  balance  technique  of  Chapter  3.  For  instance,  for 
the  tip  node  (node  1) : 


T1  -  Ti 
Rli 


+  clsTl 


(4.1) 


Here  RF1  denotes  the  thermal  resist*^®  for  convection  to 
the  fin,  whereas,  Rxi  denotes  the  thermal  resistance  for 
conduction  between  nodes.  The  letters  C  and  s  represent 
thermal  capacitance  and  the  Laplace  variable,  respectively. 
Tr2  and  TR1  are  the  recovery  temperatures. 

The  temperature  at  node  1  can  therefore  be  expressed  in 
terms  of  the  adjacent  nodal  temperatures  as  follows: 


*1 


%l 

1+r  XS 


Rpi  Ru 


(4.2) 
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where  RN1  is  given  by: 


1 


1 

%1 


rfi 


and  the  nodal  time  constant  is  given  by: 


(4.3) 


r i  —  RN1^1  (4.4) 

Subsequent  energy  balances  at  each  of  the  remaining  nodes 
yields  similar  results,  thereby  providing  sufficient 
information  to  determine  each  nodal  temperature  as  a 
function  of  time.  These  nodal  relationships  are  as  follows: 


%2  TR2  |  Ti  |  T3 

1+r  2S  Rp2  Ri2  R23  (4.5) 
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whfere,  —  =  — —  +  —  +  — -  ,  and  r  2  »  RN2C2. 
rN2  rF2  Ri2  R23 

t3  38  rN3  T2  +  Tg 

l+r  3s  r23  R3g  (4.6) 


111 

where,  ___  »  _  +  _  ,  and  r3  “  RN3C3  . 

RN3  R23  R3g 

Initially,  the  entire  system  is  assumed  to  be  at  thermal 
equilibrium  such  that  all  nodal  temperatures  are  equal  to 
the  ambient  temperature  (ground)  of  the  environment.  The 
firing  of  the  rocket  motor  disturbs  this  equilibrium  and  the 
subsequent  changes  in  the  two  recovery  temperatures  provide 
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the  inputs  to  the  resulting  heat  transfer  process.  The 
development  of  these  thermal  inputs  as  a  function  of  the 
rocket  motor  thrust  levels  is  discussed  in  the  next  chapter 
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v.  ESTIMATES  OF  HEAT  TRANSFER  PARAMETERS 


A.  GENERAL  DISCUSSION 

As  described  in  the  previous  chapter,  only  two  inputs 
are  required  to  drive  the  thermal  model  as  presently 
configured.  These  forcing  functions  are  the  result  of  the 
heating  occurring  at  the  vane  tip  and  along  the  downstream 
surface  of  the  vane.  It  should  be  noted  that  the  cooling 
effects  of  radiation  and  ablation,  both  important  in  the 
actual  situation,  are  not  yet  addressed  at  this  point  in  the 
study . 

Estimates  of  the  various  heat  transfer  parameters 
presented  in  this  chapter  are,  in  large  part,  the  result  of 
extensive  research  previously  conducted  (Ref.  8].  Most  of 
the  computations  performed  in  these  preliminary  studies  were 
derived  from  analytical  methods  typically  employed  in  the 
study  turbulent  flows  and  rocket  engine  heat  transfer 
[Refs.  9,10]. 

B.  THE  PRANDTL  NUMBER 

A  significant  feature  of  high  speed  rocket  exhaust 
flow  is  the  large  difference  in  temperature  the  gas 
experiences  in  decelerating  near  the  body  surface*  Hence, 
it  is  necessary  to  account  for  the  temperature  dependency 
of  the  gas  properties.  In  such  cases  both  viscosity  and 


2? 


thermal  conductivity  are  considered  to  vary  significantly 
with  temperature.  The  Prandtl  number  and  specific  heat, 
on  the  other  hand,  are  considered  constant  inasmuch  as 
their  variations  with  temperature  are  of  a  lower  order 
of  magnitude  than  those  of  the  other  gas  properties 
[Ref.  9:  p.  304]. 

An  important  key  to  evaluating  rocket  engine  heat 
transfer  is  the  non-dimensional  Prandtl  number  (Pr) .  The 
Prandtl  number  of  the  exhaust  gas  is  derived  from  the 
viscosity  (m  ) ,  the  thermal  conductivity  (K) ,  and  the 
specific  heat  at  constant  pressure  (Cp)  as  follows: 


K  (5.1) 

Based  on  the  above  assumption  that  Pr  and  Cp  are  constant, 

Pr  n 

—  ■  — —  -  constant 

Cp  K  (5.2) 

and  the  thermal  conductivity  must  vary  in  the  same  way  as 

the  viscosity.  Thus  an  estimation  of  the  gas  viscosity 

leads  to  a  thermal  conductivity  for  a  given  Prandtl  number 

and  specific  heat.  A  simple  equation  relating  viscosity  to 

thermal  conductivity  of  the  gas  is 

K  -  £tiCy  <5. 3) 

where  Cy  is  the  specific  heat  of  the  gas  at  constant  volume 

and  f  is  a  constant.  The  Eucken  formula  (Ref.  10 :p.  139) 

relates  f  with  the  ratio  of  specific  heats  (k)  of  the  gas, 

f  -  (9k-5)/4  (5.4) 
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from  which  it  follows  that, 


4k 

Pr  =  - 

(9JC-5)  (5.5) 

which  is  the  equation  used  in  this  work.  If  the  specific 
gas  constant  (Rg)  is  given,  the  specific  heat  at  a  constant 
pressure  is 


CP 


Rg* 

(k-l) 


(5.6) 


Also  used  in  this  study  is  the  standard  power-law  type 
formula  [Ref.  ll:p.  28]  for  determining  viscosity, 


Mref 


(5.7) 


where  n  *  0.7  (Ref.  10 :p.  8]  and  the  reference  viscosity  has 
been  taken  to  be  n  »  4.0  x  10  ~5  N-s/m2  at  Tref  ■  1000  K. 


C.  STAGNATION  POINT  HEAT  TRANSFER 

In  the  analysis  of  stagnation  point  heating,  the 
solution  to  the  boundary  layer  equation  requires  the  local 
fluid  acceleration  (B)  in  the  vicinity  of  that  point  be 
known.  With  supersonic  flow,  a  can  be  approximated  based  on 
the  assumption  that  Newtonian  flow  prevails  between  the  bow 
wave  and  the  body  [Ref.  9ip.  386].  An  approximation  for  the 
local  fluid  acceleration  has  been  derived  by  Nunn 
(Ref.  8:  p.  18]  and  is  given  byt 

fl  «  (du./dX)^  -  (U/0)[8(P./Poy)(To#/Tw)]0*5  (5.8) 
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Here  the  subscript  )  refers  to  freestream  conditions,  (Q) 
refers  to  stagnation  conditions  and  (y)  denotes  the  condi¬ 
tions  downstream  of  the  normal  shock.  The  pressure  and 
temperature  ratios  in  the  above  equation  are  known  functions 
of  the  freestream  Mach  number. 

The  Stanton  number  (St)  can  be  calculated  from  the 
following  equation  (Ref.  9:p.  366]: 

st0  *  0.57(,3D/U)0-5Pro~0’6Reo''0*5  (5.9) 

where  p  is  defined  above  and  Re  is  the  Reynolds  number. 

D.  TURBULENT  BOUNDARY  LAYER  CONVECTION 

The  analysis  of  the  heat  transfer  processes  in  high 
speed  compressible  flows  is  complicated  by  the  fact  that  a 
considerable  transfer  of  kinetic  energy  results  from  the 
deceleration  occurring  in  the  boundary  layer.  This  results 
in  temperatures  within  the  boundary  layer  that  are  in  excess 
of  the  freestream.  As  a  result  of  this  phenomenon  the  term 
" recovery  temperature11  (TR)  is  introduced  and  is  defined  as 
the  temperature  that  the  surface  will  assume  in  the  absence 
of  heat  transfer  (Ref.  10:p.  14].  Thus  the  Mrecovery 
factor*  (r)  is  defined  as  follows: 

Tr-Tw 

T0-T„  (5.10) 
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It  has  been  determined  by  Van  Driest  [Ref.  10:p.  14]  and 
others  that  for  a  turbulent  boundary,  as  is  assumed  here, 
recovery  factor  can  be  related  to  the  Prandtl  number  by 

r  =  Pr1/3  (5.11) 

A  more-difficult  problem  arises  from  taking  into  account 
the  gas  property  variation  through  the  boundary  layer  as  a 
function  of  temperature.  Fortunately,  for  compressible 
flows,  Eckert  (Ref.  10 :p.  8]  provides  a  “reference  teroper- 
ature1'  (Tref)  for  the  boundary  layer  at  which  the  various 
gas  properties  should  be  evaluated.  Eckert's  reference 
temperature  as  used  in  this  study  is  calculated  by 

Tref  *  °«5Twall  +  °*28T«  +  0*27TR  (5.12) 

Initial  calculations  have  been  based  upon  the  assumption 
that  the  wall  temperature  can  be  adequately  represented  by 
the  mean  of  the  recovery  and  ambient  temperatures. 

With  Tref  thus  defined,  the  gas  viscosity  and  thermal 
conductivity  can  be  estimated  as  described  above.  The 
Stanton  number  for  turbulent,  compressible  flow  is  then 
given  by  (Ref.  10 J: 

St  -  0.0296Pr~0*67Re*'0,2  (5.13) 

Nusselt  numbers  and  thermal  resistances  can  be  calculated 
for  the  model  nodes  affected  by  stagnation  and  turbulent 
boundary  layer  heat  transfer  from  the  Stanton  numbers 
discussed  above. 
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Thermal  resistances  are  also  dependent  upon  associated 
dimensional  parameters  (length,  area,  volume) ,  and  as  such, 
are  scale  dependent.  For  example,  the  stagnation  point 
thermal  resistance  is  proportional  to  (scale)”-1*5,  whereas 
the  turbulent  boundary  layer  thermal  resistance  decreases  as 
(scale)”1,2.  It  should  be  noted  that  the  simulations 
conducted  in  this  study  have  been  of  a  1/4 -scale  model  in 
order  to  seek  a  comparison  with  the  25%  subscale  test  data 
provided  by  NWC. 

E .  INPUT  REQUIREMENTS 

In  order  to  compute  values  for  the  thermal  resistances, 
the  stagnation  temperature,  stagnation  pressure  and  Mach 
number  of  the  freestream  flow  at  the  location  of  the  jet 
vane  must  be  determined.  Therefore,  inputs  required  for  the 
simulation  include  the  rocket  motor  chamber  pressure,  thrust 
and  characteristic  velocity  along  with  the  discharge  coef¬ 
ficient  and  pressure  ratio  of  the  nozzle.  The  propellent 
gas  properties  must  also  be  known.  These  include  the  gas 
molecular  weight  or  specific  gas  constant  and  the  ratio  of 
specific  heats. 

Given  the  inputs  mentioned  above,  the  stagnation  temper¬ 
ature  can  be  determined  by, 

T0  =  (Cd  N  c*]2/Rg  (5.14) 
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where, 

Cd  =  nozzle  discharge  coefficient 
N  «  [k(2/k+l)  (}c+1)/ (^”1)  ]0*5 
c*  =  characteristic  velocity 
Rg  =  gas  constant 
k  =  ratio  of  specific  heats 
The  Mach  number  at  the  nozzle  exit  is  computed  from 

M*2  =  [2/  (k-1)  ]  [(P^/P*)  (k'1)/k  -  1]  (5,15) 

where  Pq^/P,,  is  the  nozzle  pressure  ratio.  It  should 
noted,  that  in  this  study,  the  Mach  number  at  the  nozzle 
exit  is  assumed  to  be  that  at  the  vane. 
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VI 


A.  COMPUTATIONAL  MODEL 

As  part  of  the  preliminary  study  described  in 
Reference  5,  a  computer  program  has  been  written  which 
coded  the  3 -NODE  model  and  governing  equations  described 
in  chapter  IV,  in  the  Dynamic  Simulation  Language  (DSL) . 
Details  of  this  FORTRAN-based  language  are  included  in 
Reference  12.  The  simulation  code  is  written  such  that 
the  driving  input  is  the  rocket  motor  thrust.  The  thrust¬ 
time  profile  is  approximated  by  a  ramp  from  zero  to  a 
maximum  in  a  period  of  0.5s,  held  at  this  plateau  for 
2.5s,  and  ramped  back  to  zero  thrust  over  an  additional 
0.5s  period.  The  thrust  profile  input  used  in  the 
simulation  described  below  is  illustrated  in  Figure  6.1. 

The  magnitude  of  the  convective  heat  transfer 
resistances,  RF1  ar.d  RF2  are  a  function  of  flow  part  the 
vane  (see  Figure  4.3}.  These  resistances  are  very  large 
prior  to  ignition  and  during  burnout  and  are  minimum 
during  the  period  of  maximum  thrust.  Modeling  of  this 
behavior  is  accomplished  by  varying  the  film  coefficient 
according  to  the  thrust  schedule.  These  coefficients 
begin  and  end  at  1%  of  the  full  flow  value  corresponding 
to  the  maximum  thrust  plateau. 
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This  computational  model  incorporates  a  transfer 
function  furnished  by  NWC  to  provide  the  vane  temperature 
data  for  comparison  with  the  simulation  results  [Reference 
8].  The  transfer  function,  derived  from  experimental  test 
data,  is  given  by 


Thrust  (1.233S  +  1) (50.76S  +  1)  (6.1) 

It  should  be  noted  that  in  this  equation  the  thrust  is  in 
newtons  and  the  temperature  in  degrees  Fahrenheit  above 
ambient . 

Simulation  inputs  and  associated  heat  transfer 
parameters,  for  the  case  herein,  are  listed  in  Appendix  A 
with  the  simulation  results  presented  below. 

B.  VANE  THERMAL  RESPONSE 

DSL  simulation  results  predicting  the  vane  thermal 
response  to  input  data  of  Appendix  A.  is  shown  in  Figure 
6.2.  The  figure  shows  that  the  computed  shaft  temperature 
profile  (node  3)  is  significantly  higher  than  that 
obtained  in  the  NWC  tests. 

The  greatest  uncertainties  associated  with  the  3 -NODE 
model  result  from  the  calculation  of  "beet  estimateM 
values  for  RF1  and  Rp2  and  corresponding  film 
coefficients.  An  additional  uncertainty  arises  when 
attempting  to  adequately  model  the  mount  “heat  sink*' 


36 


(NiAiax)  aamvaaawai 


37 


Figure  6.2.  Vane  Thermal  Response 


phenomena  by  a  single  mount  resistance,  R3G.  Preliminary 
experimentation  with  simulation,  as  detailed  in  reference 
5,  has  led  to  the  insight  that  the  dominant  factor 
affecting  the  maximum  shaft  temperature  is  the  mount 
thermal  resistance.  The  dynamic  temperature  rise  of  the 
shaft,  however,  is  mainly  controlled  by  the  convective 
thermal  resistance. 

Figure  6.J  shows  the  sensitivity  of  the  shaft 
temperature  response  to  the  adjustments  of  convective  film 
coefficients.  This  figure  clearly  illustrates  that  in 
order  to  obtain  reasonable  agreement  between  simulation 
actual  it  will  be  necessary  to  lower  these  coefficients  by 
a  reduction  factor  (Rp)  of  approximately  70%.  It  has  been 
postulated  that  this  reduction  factor  may  provide  some 
insight  into  the  relative  magnitude  of  the  cooling  effects 
of  radiation  and  ablation,  which  thusfar  have  been  omitted 
rom  the  analysis.  It  is  equally  apparent  from  this 
figure,  that  the  model  is  capable  of  accurately 
reproducing  the  main  transient  characteristics  of  the  vane 
thermal  response. 

C.  SUMMARY 

The  3-NODE  model  incorporating  a  mount  thermal 
resistance  of  R3G  **  10.3  K/W  and  a  reduction  factor  of  Rp 
*  o,72  yields  results  of  the  sort  shown  in  Figure  6.4.  It 
is  noteworthy  that  the  results  of  the  simulation  indicate 
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Figure  6.3.  Shaft  Temperature  Profiles  for  Various 
Values  of  R*? 


o 


(NUiax)  aanxvaadwai 


40 


Figure  6.4.  Vane  Thersal  Response  with  RF=0.72 


I 


that  it  might  be  feasible  to  predict  the  bulk  thermal 
behavior  of  specified  critical  elements  of  the  jet  vane, 
namely  the  tip,  fin  and  shaft.  Also  important  is  the 
model  simplicity,  a  prerequisite  for  system-identification 
analysis  which  is  the  subject  of  the  next  chapter. 
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VII 


A.  FOUR  NODE  CONFIGURATION 


The  goal  of  the  work  presented  herein  has  been  to  extend 
the  results  reported  in  Reference  8  by  utilizing  additional 
experimental  data  available  from  NWC  tests.  These  include 
subscale  data  from  the  thermal  response  of  two  thermocouples 
located  in  the  shaft  and  mount  regions. 

To  fully  exploit  this  additional  information  it  has  been 
necessary  to  enhance  the  thermal  model  with  an  additional 
calculation  node.  Hence,  the  model  illustrated  in  Figure 
4.3  has  been  revised  as  shown  in  Figure  7.1. 

Energy  balance  equations  for  this  configuration  are  as 
follows: 
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Figure  7*1  Simple  Four  Node  Model  Configuration 
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B.  THE  "SYSTEM  BUILD"  MODEL 

The  system  of  equations  defined  above,  have  been  modeled 
in  the  SYSTEM  BUILD  [Ref.  13]  format  of  the  personal 
computer  version  of  the  software  package  MATRIXx  [Ref.  14]. 
The  SYSTEM  BUILD  model  is  illustrated  in  Appendix  B  and  the 
reader  should  note  the  hierachical  nesting  of  the  "super- 
blocks".  A  brief  description  of  this  arrangement  will  be 
given  here. 

To  help  illustrate  this  SYSTEM  BUILD  arrangement,  it  is 
useful  to  define  the  following  parameters: 

all  =  <a12  +  a13  *  bll) *  a12=  1/C1R12 


a21  53  l/c2*12'  a22  =  (a21  +  a23  +  ^22^  a23  =  VC2R23 
a32  ”  I/C3R23 »  333  a  (832  +  a34)^  a34  *  VC3R34 

a43  31  R4g/  (R34  +  r4g) 

bu  "  b22  -  1/C2RF2  (7.10) 

With  the  above  parameters  thus  defined,  the  state  equations 
can  be  simplified  as  follows 

T1  *  "allTl  +  a12T2  +  bllTRl  (?‘u) 

T2  “  a2lTl  “  a22T2  +  a23T3  *  b22TR2  (7.12) 

• 

T3  *  a32T2  “  a33T3  *  a34’f4  (7.13) 


In  this  study,  all  temperatures  have  been  referred  to 
ambient  temperature  (TG) ,  hence  equation  7.9  can  be 
written  as 

T4  *  a^jTj  (7.14) 
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In  this  study,  the  SYSTEM  BUILD  thermal  model  is 
comprised  of  eight  super-blocks: 

N0D1IN,  NOD2IN,  N0D3IN 
N0DE1,  NODE2,  NODE  3,  NODE 4 
VANE 

The  input  to  the  model  is  the  ramp-up,  plateau,  ramp- 
down  profile  previously  described  in  Chapter  VI.  This  input 
with  a  maximum  plateau  level  of  unity  is  supplied  to  both 
NOD1IN  and  NOD2IN,  whereas  N0D3IN  does  not  require  the  time- 
varying  input.  These  first  three  blocks  compute  the 
coefficients  necessary  for  the  state  equations  modeled  by 
the  next  three  blocks.  Super-block  N0DE4  requires  only  one 
input  from  NODE3 .  In  NODIIN,  for  example,  the  input  is 
multiplied  by  bl;L  and  TR1.  Also  in  NOD1IN  the  parameter  a12 
is  generated  by  means  of  a  step  function  and  is  summed  with 
bn  to  form  a^.  Thus  the  outputs  of  N0D1IN  are  a13_,  a12 
and  b22TR1  and  a32,  a33,  a34  respectively. 

As  illustrated  in  Appendix  B,  the  first  three  super¬ 
blocks  listed  above  are  nested  within  the  next  three  blocks 
respectively  and  NODE 4  requires  no  nested  block.  In  NODEl, 
for  example,  the  outputs  of  N0D1IN  are  combined  with  the 
external  input  T2  to  form  the  state  equation 

*1  *  ~allTl  +  a12T2  +  bllTRl  (7.15) 

* 

An  integrator  is  then  used  to  convert  to  Tj^,  which  is 
the  NODEl  output.  Again,  N0DE2,  N0DE3  and  NODE 4  are 


similar  in  function  and  provide  outputs  T2,  T3  and  T4.  The 
superblock  VANE  connects  all  four  of  these  blocks  together 
and  provides  the  simultaneous  solution  for  the  four  nodal 
temperatures. 

C.  SIMULATION  RESULTS/SYSTEM  IDENTIFICATION 

The  first  step  in  the  execution  of  the  SYSTEM  BUILD 
model  is  to  check  the  model  by  means  of  the  ANALYZE  command 
[Ref.  13].  Next,  one  of  six  integration  algorithms  may  be 
selected.  The  results  presented  herein  have  been  obtained 
with  the  variable  Kutta-Merson  method  which  is  an  explicit, 
fourth  order,  one-step  method.  This  method  employs  the 
largest  time  step  possible  while  remaining  within  error 
tolerances.  The  maximum  step  size  is  equal  to  the  time 
increment  specified  for  the  simulation.  The  final  step  is 
executing  the  SIM  command,  as  in  "Y  **  SIM  (T,  U)",  where  T 
is  the  time  vector  and  U  is  the  input,  both  of  which  must  be 
provided.  For  the  results  described  here,  T  is  a  151  x  1 
vector  created  by  dividing  15  seconds  into  0.1  second 
increments.  U  is  the  151  x  2  matrix  depicting  the  two  ramp- 
up,  unity  plateau,  ramp  down  profiles  previously  described. 

1.  NWC  Transfer  Function  Results 

Figure  7.2  shows  the  results  of  simulation  using  the 
SYSTEM  BUILD  model  to  predict  the  vane  thermal  response 
incorporating  “best  estimate1'  values  of  R34  »  3.0  K/W 
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TIME (SEC) 

Figure  7.2.  Comparison  of  Simulation  and  NWC  Transfer  Function  Results 


(a34=. 33 )  and  R4G  =  7.3  K/w  (a43=.709).  Data  from  the  NWC 
transfer  function  described  in  Chapter  VI  has  been  used  for 
comparison. 

The  model  was  set  up  for  system  identification, 
using  the  MAXLIKE  function  [Ref.  8,  15].  Figure  7.3  shows 
the  results  for  shaft  temperature  after  having  reduced  a34 
and  a43  to  the  MAXLIKE  specified  values  of  .1633  and  .550 
respectively.  These  results  are  clearly  excellent  and  give 
encouragement  that  the  Parametric  System  Identification 
(PSI)  capability  of  the  MAXLIKE  command  function  is  a 
feasible  design  tool. 

2.  NWC  1/4 -Scale  Test  Firings 

Event  (2)  and  event  (6)  are  names  given  to  separate 

1/4-scale  VLMB  rocket  motor  test  firings  conducted  by  NWC. 

The  only  planned  difference  in  the  two  tests  was  the  type  of 

propellant  burned  with  event  (2)  using  non-aluminized  HTPB 

(Hydroxyl -Terminated  Polybutadien)  and  event  (6)  burning  18% 

Aluminized  HTPB.  In  both  tests,  two  thermocouples  were 

attached  to  each  vane  in  the  shaft  and  mount  regions.  The 

thermal  response  of  these  thermocouples  were  recorded  and 

the  results  are  included  in  Appendix  C  along  with  the  rocket 

motor  ballistics  data  and  propellant  properties.  It  should 

be  noted,  however,  that  the  shaft  thermocouples  failed  1.4 

seconds  into  event  (6)  test  firing.  Hence  no  usable  shaft 

temperature  data  from  event  (6)  is  available  for  comparison 

purposes.  Event  (2)  stagnation  temperature  T0”2650K  and 
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Figure  7.3.  Comparison  of  Simulation  and  NWC  Transfer  Function  Results  After  PSI 


m 


in 


ambient  temperature  of  290K  yield  TR1=2360K  and  TR2=2260K. 
Event  (6)  stagnation  temperature  T0=3275K  and  ambient 
temperature  of  305K  yield  TR1=2970K  and  TR2=2870K. 

It  should  be  noted  that  the  mount  hardware  of  event  (6)  and 
event  (2)  each  use  different  type  fasteners. 

3*  &vent  (2)  Simulation 

Figure  7.4  is  the  result  of  a  Parametric  Simulation 
Identification  (PSI)  execution  using  initial  parameter 
values  of  a34  =  .1633  and  a43  =  .5500  taken  from  the  pre¬ 
vious  results  of  figure  7.3.  The  MAXLIKE  values  obtained 
are  a34  =  .17  and  a43  «  0.56  as  the  result  of  comparing 
model  and  experimental  shaft  temperatures.  Although  figure 
7.3  indicates  excellent  agreement  between  predicted  and 
actual  shaft  temperatures,  the  time  constant  of  NODE  4,  as 

modelled,  is  significantly  less  than  that  observed  in  actual 
test  firings. 

Eymt  ($)  Simulation  Result 

Figure  7.5  is  the  result  simulating  the  vane  thermal 
response  at  the  elevated  input  temperatures  of  TR1  »  2970K 
and  tR2  -  2870K  with  parameters  a34  and  a43  remaining  at  the 
previous  values  of  0.17  and  0.56  respectively.  This  figure 
clearly  indicates  the  disparity  between  the  model  NODE  4  and 
actual  mount  time  constants.  Figure  7.6  is  the  result  of 
increasing  r4G  in  equation  (7.9)  to  a  value  of  9.9  K/w  to 

more  closely  approximate  the  actual  mount  temperature  at 
time  *  15s. 
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Figure  7.4.  Comparison  of  Simulation  and  Event  (2)  Results 


TIME (SEC) 

Figure  7.5.  Comparison  of  Simulation  and  Event  (6)  Results,  RG=7.2  K/W 


Figure  7.6.  Comparison  of  simulation  and  Event  (6)  Results,  Rq-9.8  K/W 


From  the  results  of  Figures  7.4  and  7.5  it  became 
apparent  that  the  simple  conduction  path  linking  nodes  3  and 
4  is  an  inadequate  representation  of  the  vane/mount 
response.  It  is  obvious  that  the  T4  node  responds  with  a 
time  constant  significantly  different  than  that  of  the  T3 
node.  Correction  of  this  modeling  discrepancy  is  the 
subject  of  the  next  chapter. 
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VIII 
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V 


.  REVISED  FOUR-NODE  MODEL 

A.  REVISED  FOUR  NODE  CONFIGURATION 

In  order  to  vary  the  time  constant  at  node  4  it  was 
necessary  to  enhance  the  simple  four  node  model  with  an 
additional  capacitance  C4  located  at  this  node.  Thus,  the 
model  illustrated  in  Figure  7.1  has  been  revised  as  shown 
in  Figure  8.1. 

Energy  balance  equations  for  this  configuration  remain 
unchanged  from  those  previously  described  except  for  node 
4,  which  is  given  by 


1 

1 

1 

C4R34 

T4 

C4R34 

• 

C4R4G 

The  above  expression  can  be  simplified  as  follows 

T4  =  a43T3  “  a44T4  (8*2) 

where 

a43  "  I/C4R34 ,  a4G  ■  I/C4R iq,  a44  -  *43  +  a4g  (8.3) 

B.  REVISED  "SYSTEM  BUILD"  MODEL 

The  revised  SYSTEM  BUILD  thermal  model  is  comprised  of 
nine  super-blocks. 

NODI IN,  N0D2IN,  N0D3IN,  N0D41S 
N0DE1 ,  NODE2 ,  N0DE3 ,  N0DE4G 
VANE7 
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this  model  arrangement  is  identical  to  the  initial 

configuration  illustrated  in  Appendix  B  except  for  the 

addition  of  block  N0D4IN  shown  in  Figure  8.2  and  the 

replacement  block  N0DE4  with  block  NODE 4 G  illustrated  in 

Figure  8.3.  The  input  to  N0D4IN  is  the  external  input  T3 

from  N0DE3 .  In  N0D4IN  the  parameters  a43  and  a4G  are 

generated  by  a  step  function.  Parameter  a43  is  multiplied 

by  the  input  T3  to  form  a43T3  and  summed  with  a4G  to  form 

a44.  In  N0DE4G  the  outputs  of  N0D4IN  are  combined  with  the 

internal  input  T4  to  form  T4.  An  integrator  is  then  used 
« 

to  convert  T4  to  T4,  which  is  the  N0DE4G  output. 

Super-block  VANE7  connects  the  four  NODE  super-blocks 
together  as  illustrated  in  Figure  8.4  and  provides  the 
simultaneous  solution  for  the  four  nodal  temperatures. 

C.  SIMULATION  RESULTS/SYSTEM  IDENTIFICATION 

i.  fivont  f.<LL.BgavU-t.a 

Figurs  8.5  is  the  result  of  simulation  using 
initial  values  of  R34  *  5.9K/W,  R4G  *  7.2  K/W  and  C4  «  1.0 
J/K,  which  in  turn  yield  parameter  values  c34  -  .17, 
a43  *  *17  and  *44  “  «14»  This  shows  excellent 

agreement  between  the  model  and  actua)  mount  thermal 
responses. 

Figure  0.5  is  the  result  of  a  three-step  PS I 
execution  using  initial  parameter  values  of  a34  *  0.17, 
a43  *  0.17  and  a44  *  0.14.  In  the  first  step  only 
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Figure  8.4.  SYSTEM  BUILD  Super-Block  VANE7 
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Figure  8.5.  Comparison  of  Simulation  and  Event  (2)  Results 


«  .P"  w 


parameter  a34  is  varied  in  a  PSI  run  of  three  iterations. 
The  MAXLIKE  value  obtained  based  on  T3  is  a34  =  0.12.  In 
the  next  step,  both  parameters  a43  and  a4G  are  allowed  to 
vary  in  the  PSI  execution  on  both  nodal  temperatures  T3  and 
T4.  Although  the  MAXLIKE  function  is  capable  of  operation 
on  multiple  outputs,  this  PSI  effort  failed  due  to 
inadequate  computer  memory  storage.  As  an  alternative 
approach,  again  with  a43  and  a4G  as  the  variable  para¬ 
meters,  a  PSI  run  was  conducted  based  on  T4  alone  and 
the  MAXLIKE  values  obtained  are  a43  *  0.17  and  a4G  =  0.14. 
These  values  give  the  excellent  results  shown  in  Figure 
8.6. 

3.  SyeaL,..I3I  Jesuits 

Figure  8.7  is  the  result  simulating  vane  thermal 
response  at  the  elevated  input  temperatures  of  TR1  »  as-'OK 
and  TR2  =  2870K  with  parameters  a34,  a43  and  a4G  remaining 
constant  at  the  values  identified  in  the  previous  PSI 
analysis. 

Figure  8.8  is  the  result  of  a  PSI  execution  using 
initial  parameter  values  of  a43  *  0.17  and  a4G  *  0.14  taken 
from  the  previous  results  of  Figure  8.6.  The  MAXLIKE 
values  obtained  are  a43  *  0.1693  and  a4G  -  0.101. 

! 
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Figure  8.6.  Comparison  of  Simulation  and  Event  (2)  Results  After  PSX 
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Figure  8.7.  Comparison  of  Simulation  and  Event  (6)  Results 
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The  node  4  time  constant,  resulting  from  addition 
of  capacitance  C4  provides  a  good  representation  of  mount 
thermal  response,  as  indicated  in  the  results  outlined 
above.  These  results  are  further  discussed  in  detail  in 
the  next  chapter. 


IX.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  SUMMARY  OF  RESULTS 

From  the  results  presented  in  chapter  VIII,  it  .is 
quite  apparent  that  the  revised  four-node  model  illustrated 
in  Figure  8.1  provides  an  adequate  representation  of  the 
vane/mount  thermal  interface. 

The  values  of  parameters  a34=0.12  and  a43=0.17  obtained 
from  parametric  system  identification  analysis  using  event 
(2)  test  data  are  identical  to  the  values  resulting  from 
PSI  execution  of  event  (6)  test  data.  These  results  are 
consistent  with  the  fact  that  the  jet  vane  employed  in 
event  (2)  was  identical  to  the  vane  used  in  event  (6) . 
Similarly,  the  differences  in  MAXLIKE  resultant  values 
a4G*o.l4  and  a4G«0.l0  for  events  (2)  and  (6)  respectively 
are  an  indication  of  actual  variations  in  the  mount 
hardware  assemblies  previously  mentioned  in  Chapter  VII. 
These  results  yield  thermal  model  values  of  capacitance 
C4*o.7l  j/K  and  resistance  R34»8.3  K/W  with  R4G*10.1  K/W 
for  the  event (2)  mount  and  R4G»13.9  K/W  for  the  mount 
arrangement  used  in  event  (6) . 


B.  CONCLUSIONS 


The  goal  of  the  Modeling  and  Simulation  study  reported 
here  has  been  to  develop  a  dynamic  thermal  model  of  jet 
vane  TVC  devices  and  investigate  the  feasibility  of  using 
system  identification  methods  to  develop  tools  for  use  in 
the  design  of  such  vanes. 

The  SYSTEM  BUILD  capability  of  MATRIXx  has  been  proven 
to  be  an  excellent  computational  tool  that  provides  both 
"state-of-the-art"  graphical  model  building  and  simulation. 
This  capability  replaces  the  tedious  DSL  computer  language 
programming  employed  in  previous  modeling  and  simulation 
development.  SYSTEM  BUILD,  used  in  conjunction  with 
MATRIXx  interpreter,  allows  the  designer  to  concentrate  his 
or  her  efforts  on  design,  analysis  and  simulation  in  an 
efficient  graphical  environment. 

The  development  of  MATRIXx  as  a  design  tool  is  further 
manifested  by  its  system  identification  capabilities,  in 
particular  the  MAXLIKE  function  used  in  the  work  reported 
here.  In  this  study,  the  power  of  parametric  system 
identification  has  been  exploited  in  removing  the  some  of 
the  uncertainties  associated  with  the  revised  four-node 
model . 
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The  results  presented  in  this  thesis  yield  a  measure  of 
confidence  that  the  thermal  behavior  of  the  jet  vane  TVC 
device  can  be  modeled  using  a  relatively  simple  "lumped- 
capacitance"  type  model  configuration. 

C.  RECOMMENDATIONS 

Further  simulation  and  PSI  analysis  should  be  performed 
using  a  computer  system  such  as  a  VAX-2000  work  station  and 
thereby  ensuring  adequate  storage  capacity  is  available  to 
fully  exploit  the  multi-input/multi-output  systems 
capability  of  the  MAXLIKE  algorithm. 

Also,  a  FORTRAN  compiler  should  be  incorporated  into 
whatever  system  MATRIXx  is  loaded,  in  order  to  enhance  the 
versatility  of  inputting  test  data.  Data  generated  by 
other  programs,  or  obtained  from  experiments  must  be 
converted  to  or  stored  in  a  form  which  can  be  read  by  the 
MATRIXx  LOAD  command .  MATSAV  is  a  FORTRAN  subroutine 
designed  to  save  external  data  in  a  proper  format. 

A  logical  next  step  in  this  study  would  be  to  adjust 
the  parameters  of  the  present  model  to  full-scale  values 
using  direct  scaling  procedures  such  as  those  given  in 
chapter  V.  A  comparison  between  the  simulated  thermal 
response  of  the  model  and  recently  acquired  full-scale  test 
data  could  be  made  and  additional  PSI  conducted  if 


necessary 


APPENDIX  A 


PRELIMINARY  DSL  SIMULATION  PARAMETERS 
ROCKET  MOTOR 

Cd  =  Nozzle  Discharge  Coefficient  =  .934 
r  =  0.65 

C*  =  Characteristic  Velocity  =  1512  m/s 
Rg  =  Gas  Constant  =  318.5  ra2/s2  “  K 
K  =  Ratio  of  Specific  Heats  =  1.21 
T0  =  Stagnation  Temperature  =  2650K 
P0®  =  Maximum  Chamber  Pressure  =  15.76MPa 
Pqco/Pco  a  Nozzle  Fressure  Ration  =  186 

M«  =  Mach  Number  at  Nozzle  Exit  =3,75 
Maximum-  2325N 

SXASB&XIQN  POINT 

Staton  Number  »  8.53  x  10~3 
Nusselt  Number  -  43.1 
Film  Coefficient  -  6.46  x  103  W/M2  -  K 
Thermal  Resistance  -  5.69  K/w 

Stanton  Number  *  3.01  x  10~3 
Nusselt  Number  *  438 
Film  Coefficient  -  2.25  x  103  W/M2  -  K 
Thermal  Resistance  *  0.634  K/w 


APPENDIX  C 


NWC  TEST  FIRING  PARAMETERS 


Event ( 2 ) 

Event  ( 6 ) 

Propellant 

0%  AL  HTPB 

18%  AL  HTPB 

Molecular  Weight 

26.1 

30.95 

K  = 

1.210 

1.145 

C*  - 

4930  ft/sec 

5175  ft/sec 

II 

O 

.934 

.934 

Burn  rate  Coefficient  = 

.05703  in/sec 

.00807  in/sec 

Burn  rate  exponent  = 

.28063 

.54777 

Ambient  Conditions 

70°  F,  13.7  PSIA  70° 

F,  13.7  PSIA 

Propellant  Mass  = 

6.17  lbm 

6.63  lbm 

Propellant  Density  = 

.06185  lbm/ in3 

.06671  lbm/in 

Initial  Throat  Area  * 

.14750  in2 

.26210  in2 

Post-fire  Throat  Area* 

.15135  in2 

.26558  in2 

Max  Pressure  = 

2286.8  PSI 

1078.5  PSI 

Max  Thrust  = 

514  Ibf 

420  lbf 

Total  Impulse  » 

1535  lb-sec 

1667  lb-sec 

ISP  = 

248.8  sec 

251.4  sec 

Exit  Dia  <* 

1.8125” 

1.815” 

Max  Exit  Pressure  • 

12.3  PSI 

14.8  PSI 

8C 
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